How well do we know the neutron structure function? 
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We present a detailed analysis of the uncertainty in the neutron F2„ structure function extracted 
from inclusive deuteron and proton deep-inelastic scattering data. The analysis includes experi- 
mental uncertainties as well as uncertainties associated with the deuteron wave function, nuclear 
smearing, and nucleon off-shell corrections. Consistently accounting for the Q 2 dependence of the 
data and calculations, and restricting the nuclear corrections to microscopic models of the deuteron, 
we find significantly smaller uncertainty in the extracted F-i n fF^ v ratio than in previous analyses. 
In addition to yielding an improved extraction of the neutron structure function, this analysis also 
provides an important baseline that will allow future, model-independent extractions of neutron 
structure to be used to examine nuclear medium effects in the the deuteron. 



Because the free neutron is an experimentally im- 
practical scattering target, extracting its structure func- 
tion, F 2n , requires cross section data from inclusive 
deep-inelastic scattering (DIS) measurements on proton 
and deuteron targets, together with a model describing 
the smearing produced by the nuclear binding in the 
deuteron. Previous extractions [1-3 using a variety of 
models of the deuteron bound state yielded a large range 
of FmjFip values from the same proton and deuteron 
data, indicating large theoretical uncertainties, particu- 
larly at high values of the Bjorken variable x — see Fig.[T] 
In addition to limiting our ability to extract the neu- 
tron structure function, the large spread of results, even 
among extractions including only traditional nuclear ef- 
fects such as Fermi motion and binding, has made it dif- 
ficult to identify a reliable baseline which could be used 
to search for more "exotic" nuclear effects such as the 
modification of the nucleon structure function in nuclei 
or non-nucleonic degrees of freedom. 

A more recent extraction of the neutron i^n structure 
function [3] showed that some of the variation in results 
can be attributed to inconsistent treatment of kinematics 
of the data and calculations. In particular, it was vital 
to properly account for the Q 2 dependence of the proton 
and deuteron data, especially for x > 0.7 [IJ[5]. 

Accurate information on the FmjFi v ratio at large 
x is essential for a number of reasons. These include 
constraining leading-twist parton distribution functions 
(PDFs) in the region x > 0.7, which are an important 
input for QCD background calculations in searches of 
physics beyond the standard model at the Tevatron and 
the LHC, and in neutrino oscillation experiments. The 
ratio F 2n /F2 P also provides insight into the nonpertur- 
bative quark-gluon dynamics in the nucleon [21 [5] . There 
are several predictions for F^nJF^p based on symmetry 
arguments that determine the dominant contributions as 
x — > 1, ranging from F^n/Fap — 2/3 for exact spin-flavor 
SU(6) symmetry, to 3/7 assuming helicity conservation 
through hard gluon exchange [7,, to 1/4 when SU(6) sym- 
metry is broken through scalar diquark dominance [HI [5] ■ 
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FIG. 1: (Color online) Previous extractions of Fm/F2 P , 
using microscopic deuteron calculations (filled symbols) 
or extrapolations of nuclear effects in heavier nuclei (open 
symbols): (black) triangles from Ref. pQ, (blue) squares from 
Ref. [2j, (red) circles from Ref. [3J. The arrows indicate 
theoretical predictions for the x — > 1 limit (see text). 



In this paper we extend the analysis of Ref. [I] by per- 
forming a detailed study of the model dependence of the 
extraction procedure, systematically assessing the uncer- 
tainties arising from the deuteron wave function at short 
distances, nucleon off-shell effects, and different nuclear 
smearing models used to compute the nuclear corrections. 
The goal is to determine the degree to which the F^n neu- 
tron structure function can be determined at large x, us- 
ing a consistent treatment of input data sets and realistic 
models of nuclear effects, and employing a methodology 
that is transparent and conceptually accessible. 

Following Ref. [4], we consider F2d/F 2p data from 
SLAC, BCDMS, and NMC, as compiled in Ref. [1] for 
3 < Q 2 < 230 GeV 2 , using the full Q 2 range to de- 
termine the Q 2 dependence of the ratio. The extrac- 
tion of F2n/F2 P itself is limited to data in the range 
8 < Q 2 < 32 GeV 2 , and the results are interpolated 
to a fixed Q 2 = Q\ = 16 GeV 2 . The interpolation to 
Qq differs from the simple average over the limited Q 2 
range by at most 0.7%, with a typical correction of 0.3%. 
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Most previous extractions of F 2n /F 2p used of F 2 d/F 2p 
from the analysis of Ref. [I], in which Q 2 varies from 4.7 
to 23.6 GeV 2 over the x range of the data. However, the 
extractions treat the F 2 d/F 2p ratios as though they were 
all at some average Q 2 value. Such extractions neglect 
the Q 2 dependence of the nuclear effects, which have been 
found to be significant at large x [H E] . By interpolating 
all of the source data to a common Qq and extracting 
F 2n /F 2p using different models evaluated at the same 
scale, a more systematic and meaningful assessment of 
the model dependence can be made. 

The extraction of F 2n / ' F 2p proceeds by assuming that 
the deuteron structure function F 2 d can be expressed as a 
sum of smeared proton and neutron structure functions, 
and an additional correction, 5F 2 d, that goes beyond the 
convolution approximation, F 2 d = F 2p + F 2n + SF 2 d, 
where F 2 n = Sn F 2 n is the smeared nucleon structure 
function (N = p,n), and Sn is the smearing ratio. The 
term 5F 2c [ includes any corrections, such as relativistic or 
nucleon off-shell corrections [TQ1 E] i that cannot be ex- 
pressed as a convolution of a smearing function and the 
free nucleon structure function. Parametrizing this cor- 
rection as a ratio to the total deuteron structure function, 
A = SF 2c i/F 2( i, the F 2n /F 2p ratio can then be extracted 
using a modified smearing factor Sn = <SW(1 — A), 
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The neutron to proton ratio can thus be extracted from 
F 2 d/F 2p and F 2p data, and a model of the smearing ratios 
Sn and the off-shell correction A. 

To standardize comparisons of the different calcula- 
tions, all extractions use the same values for F 2p and F 2n 
to compute the smearing functions Sn(x). We use the 
parametrization of the world's F 2p data and the extracted 
neutron structure function from Ref. [4]. Each calcula- 
tion yields a slightly different F 2n , however, the impact 
of using this modified F 2n value to update the calcula- 
tion of S n is small compared to the other uncertainties 
(F 2n /F 2p changes by less than 0.01 at x — 0.85) [T2"] . 

The assessment of the model dependence of the ex- 
tracted F 2n /F 2p ratio ultimately depends on the choice 
of nuclear models used in the analysis. The introduction 
of some degree of bias is therefore inevitable, although 
we aim to make the selection criteria as objective as pos- 
sible by restricting ourselves to microscopic calculations 
involving high-precision TV TV potentials that give realistic 
descriptions of the deuteron bound state. Common fea- 
tures of the models surveyed in this analysis include the 
use of realistic deuteron wave functions which account 
for nuclear Fermi motion and binding, an exact treat- 
ment of finite-Q 2 kinematics, and allowance for possible 
nucleon off-shell corrections in the deuteron [31 01 [T31 [T3] . 
We exclude models that involve extrapolations of nuclear 
medium modifications observed in structure functions of 



heavy nuclei to the deuteron [TJ [31 [T5] , which typically 
invoke a very large nuclear density, contain no Q 2 de- 
pendence, and effectively assume S n = S p . Note that 
this has less impact than one would expect from Fig. [T] 
If one repeats these analyses taking into account the Q 2 
dependence of the F 2 d/F 2p data, the ratios are below 0.5 
at x > 0.6 for all of the extractions. 
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FIG. 2: (Color online) Neutron to proton structure function 
ratio F^njF-z-p calculated using various deuteron wave func- 
tions (see text), within the WBA smearing model [11H13] . 

The dependence on the choice of deuteron wave func- 
tions is illustrated in Fig. [2] which shows the re- 
sults for F 2n /F 2p using different modern nonrelativis- 
tic (N3LO [IB], CD-Bonn [T7], AV18 [TJ]) and relativis- 
tic (WJC-1, WJC-2 [TO]) deuteron wave functions. The 
smearing factors for each of the calculations were com- 
puted using the weak binding approximation (WBA) 
model |11[ 113] . which is derived by expanding the nu- 
cleon correlation function in the nucleus in powers of the 
nucleon momentum p up to order p 2 /M 2 , where M is the 
nucleon mass. The gray band represents the RMS spread 
of the F 2n /F 2p ratios for the five wave functions consid- 
ered, which is the la band for the wave function uncer- 
tainty if we treat each of the iVTV potentials on an equal 
footing. The ratio F 2n /F 2p clearly becomes increasingly 
sensitive to the choice of deuteron wave function at larger 
x values, reflecting the larger uncertainty in the NN in- 
teraction at short distances or high nucleon momenta in 
the deuteron. 

The dependence of F 2n /F 2p on the model used for the 
smearing function (or smearing factor) and off-shell pre- 
scription is illustrated in Fig. [3| The curves show the 
results for the given calculation, averaged over the po- 
tentials shown in Fig. [2] For calculations where not all 
potentials were available, the WBA result was used to 
extrapolate to the average potential. Note that some 
smearing calculations, such as those used in Ref. [TJ, have 
been omitted as they represent calculations similar to 
those included here, but with additional numerical ap- 
proximations. The solid curves in Fig. [3] are the results 



3 



0.6 
0.5 
a 0.4 
U~ 0.3 
0.2 
0.1 
0.0 




WBA + KP 1 .8% 
WBA + KP 1 .5% 
MST on-shell 
WBA + MST off-shell 
WBA on-shell 
ACHL on-shell 
RT on-shell 



0.4 



0.6 

x 



0.8 



1.0 



FIG. 3: (Color Online) The ratio J^n/i^p calculated using 
different smearing models, taking the average of the NN po- 
tentials from Fig. [2] The solid curves are on-shell calcula- 
tions: WBA [niCES], relativistic MST [3J [TU] , the light-front 
ACHL g], and the RT [14] models. The dashed curves show 
the WBA result with various off-shell prescriptions: MST [10] , 
and KP [11] with two different nucleon swelling parameters, 
1.5% and 1.8% (see Ref. [20]). The ordering of the curves in 
the legend is based on the value of the extraction at x — 0.9. 



of smearing calculations with on-shell nucleon structure 
functions, with their spread indicating the uncertainty 
associated with the smearing function. The WBA cal- 
culation (solid black curve) is a modern calculation that 
makes minimal approximations, and is used as the base- 
line for showing the results of calculations including nu- 
cleon off-shell corrections. The four WBA results (solid 
black curve and the three dashed curves) indicate the 
model dependence in the off-shell prescriptions. 

To estimate the combined uncertainty, we take the 
RMS spread of all of the extractions of F 2n /F 2p shown in 
Fig. [3] at each x value, indicated by the light gray band. 
Similar uncertainties are obtained if the model depen- 
dence of the smearing function and off-shell contributions 
are extracted separately and combined in quadrature. 
The uncertainty associated with the smearing function 
is essentially negligible up to x = 0.6, but is comparable 
to the off-shell corrections for x > 0.75. 

Figure [4] shows the combined uncertainty range (gray 
band) compared to the range of results shown in Fig. [I] 
(red hatched region). The central result is taken as the 
global average of the extracted F 2n /F 2p values obtained 
in Fig. [3j The individual uncertainties associated with 



the experimental systematic uncertainties (evaluated in 
Ref. (4]), the dependence on the deuteron wave func- 
tion, and the dependence on the smearing function and 
off-shell effects (labeled "Model Uncertainty" in Fig. [4]) , 
are shown separately, as well as the sum of uncertainties 
added in quadrature. 
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FIG. 4: Fzn/Fip ratio together with the individual contribu- 
tions to the systematic uncertainty and the quadrature sum. 
The red hatched region corresponds to the uncertainty range 
in Fig. |T| 



Our analysis provides a significantly narrower range 
of results than that evident in Fig. [T] for the full spec- 
trum of models. At x = 0.85, for example, the range in 
Fig. [I] spans 0.2 < F 2n jF<2,p < 0.7, whereas the present 
analysis suggests a ler range of 0.18 < F 2n /F 2p < 0.32. 
The tighter bounds are largely due to the exclusion of 
models involving extrapolation of nuclear medium ef- 
fects from heavy nuclei, and would appear to exclude the 
SU(6) predictions of F 2n /F 2p — > 2/3, while favoring the 
lower estimates consistent with the partonic lower limit 
of F 2n /F 2p -> 1/4. 

The dependence of the extraction on the choice of nu- 
clear model is further illustrated in Fig. [5j which shows 
the total uncertainty bands corresponding to two differ- 
ent subsets of results. The solid lines show the range ob- 
tained by taking only the on-shell extractions in Fig. [3j 
which yield 0.16 < F 2n /F 2p < 0.28 at x = 0.85. While 
much of the range at large x is below F 2n /F 2p = 0.25, 
which in the parton model is forbidden by the require- 
ment that PDFs are positive, the results are not incon- 
sistent with F 2n /F 2p > 0.25. This range can be thought 
of as a baseline for effects beyond the on-shell convolu- 
tion approximation, and comparison of these results to 
model-independent extractions of F 2n /F 2p can be used 
to isolate off-shell contributions or other more exotic nu- 
clear effects. 

On the other hand, most modern quantitative analyses 
of nuclear structure functions and the nuclear EMC effect 
require the inclusion of some modification of the nucleon 
structure function in a nuclear medium |111 121H23| . Re- 
stricting the set to only models that incorporate off-shell 
effects, one obtains the dashed blue band in Fig. [5j As 
expected, this leads to a higher range for the neutron 
to proton ratio, 0.25 < F 2n /F 2p < 0.36 at x — 0.85. 
Note that in both cases, the experimental systematics 
and deuteron wave function dependence is included in 
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the bands. 
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FIG. 5: F'2n/F2 P ranges for on-shell (solid) and off-shell 
(dashed) extractions. 

At leading order in the strong coupling constant, the 
nucleon structure functions are given by the charge- 
squared weighted sum of the u and d quark distributions. 
In this approximation the extracted value of F2 n /F 2p is 
directly related to the ratio of d to u quark distributions, 

d = AF 2n /F 2p - 1 

u 4 - F 2n /F 2p ' 1 ' 

The resulting d/u ratio is shown in Fig [6j along with 
the fractional uncertainty (inset), for the full range of 
models, as well as for the on-shell and off-shell models 
from Fig. [5] Such an extraction neglects higher-order 
perturbative QCD corrections, target mass corrections, 
and higher twist effects. On the other hand, higher twists 
will cancel in this ratio, unless they differ for the proton 
and neutron [21] , as will target mass corrections to a large 
extent, although some residual prescription dependence 
survives at large values of x |25j . Therefore, even though 
Eq. ^ is approximate, the results in Fig. [6] do serve to 
illustrate how uncertainties in F 2n /F 2p propagate to the 
d/u ratio. 

The absolute uncertainty on the d quark distribution 
is small compared to the overall size of the u quark dis- 
tribution (< 10% for all x values shown). In contrast, 
the fractional uncertainty on the d quark PDF is large, 
and will yield significant uncertainties on quantities sen- 
sitive to small relative uncertainties in d, such as the 
PDF inputs for cross section calculations in high-energy 
collisions. Although the error band in Fig. [6] includes a 
significant region with d < for x > 0.8, corresponding 
to F 2n /F 2p < 0.25, the la bands are never inconsistent 
with a positive PDF. 

While the extracted d/u ratio in Fig. [6] is illustra- 
tive of the reduced uncertainty from the restricted range 
of nuclear corrections considered here, the primary goal 
of the present work is an extraction of the total F 2n 
rather than a separation of the PDFs from target mass 
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FIG. 6: Ratio d/u extracted from F 2n /F 2p at Q 2 = 16 GeV 2 
using Eq. Q, with the gray band indicating the total un- 
certainties as shown in Fig. [4] with the dashed (solid) curves 
representing the on-shell (off-shell) nuclear model results. The 
inset shows the growing fractional uncertainty on d/u as the 
ratio becomes smaller for x — > 1. 

and higher twist contributions. A complementary study 
which quantified the effects of nuclear corrections on 
PDFs within a global QCD analysis was recently per- 
formed by the CJ Collaboration [20]. The largest im- 
pact of the nuclear effects, which were computed using 
the WBA smearing function with nucleon off-shell cor- 
rections from the modified KP model (see Fig. [3]), was 
found for the d quark PDF at large x. The uncertainties 
in the resulting d/u ratio were similar to those in Fig. |6j 
although with a somewhat larger spread. In particular, 
the range in Ref. [20] for the neutron to proton ratio at 
x = 0.85 was found to be 0.32 < F 2n /F 2p < 0.50 at 
Q 2 = 16 GeV 2 , where the upper and lower limits were 
obtained using models with the minimum and maximum 
nuclear corrections, respectively. This is a larger range 
than found here, although it represents the full range of 
results, rather than an estimated la error band. Tak- 
ing a linear sum of all uncertainties from Fig. [4] yields 
the range 0.04 < F 2n /F 2p < 0.42. The upper limit is in 
reasonable agreement with the CJ global fit, while the 
lower limit is significantly smaller. This is because PDFs 
in global QCD analyses are constrained to be positive a 
priori, forcing the overall d/u bands (and consequently 
F 2n /F 2p ) to lie higher than those in Fig. [6] 

In summary, we have extracted the F 2n /F 2p structure 
function ratio from a range of models of the deuteron 
structure, and estimated the uncertainties associated 
with the choice of nuclear smearing model, deuteron wave 
function, nucleon off-shell corrections, and experimental 
uncertainties. Restricting the analysis to models based 
on microscopic calculations of deuteron structure, and 
excluding those that rely on extrapolations from heav- 
ier nuclei, we find a range of results that is significantly 
smaller than in some previous extractions. The general 
consistency between different microscopic models eval- 
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uated here provides a reliable baseline for the neutron 
structure assuming only traditional nuclear effects, so 
that future model-independent extractions of the neutron 
structure function 26, 27] can both improve our knowl- 
edge of the neutron structure and have the potential to 
provide signatures for more exotic nuclear effects. 
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